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In  this  work  we  investigate  the  effects  of  hydrogen  addition  on  the  flame  structure  of  MILD 
combustion  both  experimentally  and  numerically  using  a  laminar-jet-in-hot-coflow  (LJHC) 
geometry.  The  addition  of  hydrogen  appreciably  decreases  the  flame  height  ( —  25%), 
however  only  modestly  affects  the  maximal  flame  temperature  and  the  thickness  of 
combustion  zone.  The  NO  distribution  is  dominated  by  mixing  of  the  NO  formed  in  the 
coflow  with  the  reaction  products  of  the  diluted  fuel,  with  negligible  NO  formation  from  the 
fuel  in  all  flames  studied.  The  numerical  data  are  in  reasonably  good  agreement  with  the 
measurements. 

Copyright  ©  2013,  Hydrogen  Energy  Publications,  LLC.  Published  by  Elsevier  Ltd.  All  rights 

reserved. 


1.  Introduction 

MILD  (moderate  and  intense  low-oxygen-dilution)  combus¬ 
tion  is  a  new  technology  combining  high  efficiency  and  low 
pollutant  emissions  in  industrial  heating  processes  [1],  This 
method  adopts  highly  preheated  and  diluted  air  and/or  fuel  to 
homogenize  the  temperature  field  and  minimize  NOx  emis¬ 
sion  during  combustion  processes.  An  important  aspect 
essential  for  successful  implementation  of  this  combustion 


technology  is  the  robustness  of  a  practical  system  to  varia¬ 
tions  in  fuel  composition;  under  industrial  conditions,  fuels 
may  vary  from  natural  gas  to  H2/CO  mixtures  containing  un¬ 
saturated  hydrocarbons,  and  the  combustion  systems  must 
continue  to  operate  within  the  desired  performance  envelope. 
Currently,  there  are  relatively  few  studies  considering  varia¬ 
tions  in  fuel  composition  on  MILD  combustion.  For  example,  it 
was  shown  in  [2]  that  adding  ethane  or  N2  to  methane  both 
widened  the  operating  range  for  MILD  combustion  as 
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compared  to  methane  itself.  Another  study  [3]  found  that  the 
addition  of  N2  or  C02  to  the  fuel  (a  variety  of  fuel  mixtures 
using  methane,  ethylene,  and  propane  were  investigated) 
reduced  the  NOx  emissions  and  facilitated  flameless  oxida¬ 
tion.  In  [4]  it  was  reported  that  NOx  emissions  are  not  sensitive 
to  the  hydrogen  concentration  under  MILD  conditions.  A 
higher  fuel  jet  velocity  was  required  to  achieve  the  transition 
from  the  conventional  flame  structure  to  MILD  combustion 
when  using  a  hydrogen-methane  fuel  [5]  compared  to  pure 
methane,  but  a  lower  furnace  temperature  could  be  used. 
Only  few  studies  considered  the  physical/chemical  origins  of 
the  effects  of  fuel  composition,  in  particular  H2  addition,  on 
either  the  stability  of  the  MILD  regime  or  NO  formation.  For 
example,  it  was  shown  both  computationally  and  experi¬ 
mentally  in  [6]  that  the  reaction  zone  structure  is  very  similar 
for  the  different  fuels  considered  when  hydrogen  is  added  to 
the  fuel  stream.  Hydrogen  was  found  necessary  for  the 
experimental  flames  to  stabilize.  A  numerical  study  [7] 
demonstrated  that  the  hydrogen  addition  to  methane  leads 
to  improved  mixing,  an  increase  in  the  turbulent  kinetic  en¬ 
ergy  decay  along  the  flame  axis  and  flame  entrainment, 
higher  reaction  intensities,  improved  mixture  ignitability  and 
enhanced  rate  of  heat  release.  Numerical  and  experimental 
investigations  of  burners  operating  in  MILD  combustion 
regime  and  fed  with  methane  and  methane-hydrogen  mix¬ 
tures  showed  the  need  of  a  detailed  chemistry  approach  [8] 
and  the  need  of  a  proper  turbulence/chemistry  interaction 
treatment  [9]  to  capture  the  volumetric  features  of  MILD 
combustion.  The  effect  of  the  composition  of  methane/ 
hydrogen  mixture  on  mild  flameless  combustion  was  studied 
in  [10]  on  a  laboratory-scale  facility.  In  [11]  it  was  shown  that 
hydrogen  oxidation  kinetic  significantly  interacts  with 
methane  Mild  Combustion  process.  See  also  [12]  for  more 
references  on  the  flameless  oxidation  of  various  fuels.  We 
note  that  all  above  mentioned  work  was  done  in  turbulent 
flames.  To  our  knowledge  there  are  no  combined  experi¬ 
mental  and  numerical  studies  of  the  effects  of  fuel  composi¬ 
tion  on  MILD  flames  in  good  defined  conditions  where 
chemical  models  can  be  tested. 

In  this  work  we  investigate  the  effects  of  hydrogen  addition 
on  both  spatial  structure  and  the  NO  formation  of  a  MILD 
flame  burning  using  an  axisymmetric  laminar-jet-in-hot- 
coflow  (LJHC)  burner  both  experimentally  and  numerically. 
Towards  this  end,  we  expand  our  previous  studies  performed 
on  the  MILD  flame  with  methane  as  a  basis  fuel  (this  base 
flame  was  investigated  by  the  authors  in  [13—15])  by  adding 


hydrogen  to  the  fuel.  The  distributions  of  temperature  and 
fractions  of  major  species  and  NO  are  measured  using  spon¬ 
taneous  Raman  scattering  and  laser-induced  fluorescence 
(LIF). 


2.  Experimental 

A  schematic  of  the  LJHC  burner  is  shown  in  Fig.  1.  An  upright 
ceramic  fuel  tube  (i.d.  of  7  mm,  wall  thickness  of  1.5  mm)  was 
surrounded  by  a  coflow  annulus  (a  perforated  ceramic  tile 
with  an  overall  diameter  of  58  mm).  The  experimental  con¬ 
ditions  (the  exit  velocity,  temperature  and  composition)  for 
the  fuel  and  oxidizer  are  reported  in  Table  1  for  two  Cases 
studied.  The  hot  oxidizer  stream  was  formed  by  the  com¬ 
bustion  products  of  a  nearly  one-dimensional  flame  of  a  fuel- 
lean  N2/CH4/02  mixture,  see  [14,15]  for  more  details  on  the 
operation  of  the  burner.  The  experimental  and  numerical  re¬ 
sults  for  Case  M  were  reported  and  thoroughly  described  in 
our  previous  works  [13-15],  Here  we  study  the  effects  of  the  H2 
addition  to  the  fuel  on  the  spatial  structure  of  the  MILD  flame 
(Case  M).  Methane  and  hydrogen  were  mixed  while  keeping 
the  total  flow  rate  of  the  fuel  constant.  The  flows  of  all  gases 
were  measured  using  calibrated  mass  flow  meters  (Bronk- 
horst).  The  velocities  reported  in  Table  1  are  the  flows  of 
coflow  and  fuel  (measured  at  standard  conditions)  normalized 
by  the  corresponding  exit  areas.  The  profiles  of  temperature 
and  species  fractions  were  obtained  by  moving  the  burner 
axially  and  radially  (vertically  and  horizontally,  respectively) 
by  a  precision  positioner  (Parker,  positioning  uncertainty  less 
than  0.1  mm). 

The  flame  temperature,  major  species  (CO,  C02,  N2,  H2, 
H20,  CH4  and  02)  and  NO  fractions  were  measured  using 
spontaneous  Raman  scattering  and  Laser  induced  Fluores¬ 
cence,  respectively.  The  optical  scheme  for  the  spontaneous 
Raman  measurements  is  essentially  identical  to  that 
described  in  [16].  At  the  exit  plane  of  an  f/4  Acton  Research 
Spectra  Pro  2300i  spectrometer,  a  Pi-Max  intensified 
1024  x  1024  pixel  CCD  camera  (Princeton  Instruments,  13  pm 
pixel  size)  was  mounted.  In  the  present  study,  40  pixels  were 
binned  along  the  axis  of  the  laser  beam,  integrating  the  signal 
over  distances  of  roughly  1  mm,  yielding  the  radial  spatial 
resolution  of  the  experiment.  The  temperature  measurements 
were  performed  using  moderately  resolved  Raman  spectra  of 
N2,  the  concentration  measurements  using  low-resolution 
Raman  spectra  which  included  the  spectral  features  of  all 


Table  1 

-  The  experimental  conditions. 

Case 

Coflow 

Fuel 

ua  (cm/s) 

o2b 

n2 

H20 

C02 

NO  (ppm) 

Tc  (K) 

va  (cm/s) 

ch4 

h2 

n2 

Td  (K) 

M 

15.5 

0.036 

0.732 

0.145 

0.087 

10 

1530 

5.8 

0.18 

0.82 

1150 

H 

15.5 

0.036 

0.732 

0.145 

0.087 

10 

1530 

5.8 

0.09 

0.09 

0.82 

1185 

a  The  exit  velocities  are  reported  at  Standard  conditions  (T  =  273  K,  P  =  1  atm), 
b  The  composition  of  gases  is  given  by  volume  %. 
c  Coflow  temperature  in  the  exit  plane  of  the  fuel  tube, 
d  Fuel  temperature  at  3  mm  above  the  exit  plane  of  the  fuel  tube. 
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Fig.  1  -  Schematic  of  the  laminar-jet-in-hot-coflow  burner. 


major  species  in  our  flames.  The  optical  scheme  of  the  LIF 
measurements  and  the  method  to  derive  the  absolute  NO 
concentration  are  described  elsewhere  [14]  and  in  its  refer¬ 
ences.  In  brief,  a  Sirah  PrecisionScan  tunable  dye  laser  pum¬ 
ped  by  the  second  harmonic  of  a  Spectra-Physics  Quantra  Ray 
Pro  250-10  Nd.'YAG  laser  was  used  to  generate  wavelengths 
near  226  nm,  with  pulse  duration  t  ~  10  ns,  by  mixing  the  dye 
output  with  the  third  harmonic  of  the  same  laser.  The  laser 
beam  was  focused  to  the  flame  by  a  quartz  lens  (/=  85  cm). 
The  laser  pulse  energies  were  typically  less  than  10  jiJ  to 
assure  linearity  of  the  LIF  signals.  Fluorescence  from  the 
centre  of  the  flame  was  collected  at  right  angles  by  a  quartz 
lens  (Nikkor  f/4.5)  and  focused  onto  an  Acton  Research  Cor¬ 
poration  SpectraPro  2150i  spectrometer  (f/4,  10  nm/mm).  The 
spectrometer  bandpass  was  set  to  30  nm  to  capture  the  entire 
NO  A-X(0,  3)  band.  The  fluorescence  signals  were  detected  by 
an  Electron  Tubes  9659QB  photomultiplier.  The  spatial  reso¬ 
lution  of  the  NO  measurements  along  the  laser  beam  in  this 
work  was  ~1.5  mm. 


3.  Computations 

The  mathematical  formulation  and  numerical  methodology 
used  to  perform  the  modeling  was  extensively  discussed  in 
[15],  In  short,  the  governing  equations  for  mass,  momentum, 
energy  and  species  are  written  in  the  cylindrical  coordinate 
system.  The  equations  are  solved  using  a  finite-volume,  den¬ 
sity  based,  body-fitted  multi-block  quadrilateral  mesh 


together  with  a  parallel  adaptive  mesh  refinement  (AMR) 
following  the  approach  developed  in  [17].  The  Soret  effect  was 
not  included  in  the  calculations,  see  below.  The  numerical 
calculations  were  performed  with  the  GRI-Mech  3.0  detailed 
chemical  reaction  mechanism  [18].  The  GRI-Mech  3.0  mech¬ 
anism  was  adopted  with  some  modification  of  the  NO  chem¬ 
istry,  see  [19,20].  Namely,  the  reaction  CH  +  N2  ->  HCN  +  N  in 
GRI  3.0  mechanism,  is  changed  to 

CH  +  N2  — >  NCN  +  H,  (Rl) 

while  keeping  the  same  rate  constant,  and  4  reactions 
describing  the  reactions  of  NCN  with  O,  H,  OH,  02  are  added: 


NCN  +  O  — >  CN  +  NO  (R2) 

NCN  +  H  ->  HCN  +  N  (R3) 

NCN  +  OH  — ►  HCN  +  NO  (R4) 

NCN  +  02  ->  NCO  +  NO  (R5) 


The  calculations  were  performed  using  rate  coefficients  for 
these  four  reactions  taken  from  Lin  and  co-workers  [21-24]; 
the  values  used  are  summarized  in  [19].  Transport  is  described 
by  a  mixture  averaged  model  in  which  the  diffusion  velocity  of 
each  individual  species  is  calculated  assuming  Fickian 
diffusion. 

Here  we  investigate  the  importance  of  the  Soret  effect  in 
the  computations  of  Case  H  by  comparing  the  chemical  time 
scales  of  corresponding  ID  igniting  counterflow  diffusion 
flame  (see  Fig.  2).  These  ID  flames  are  computed  with  the 
same  composition  and  temperature  of  fuel  and  oxidizer  of  the 
Case  H2  using  the  GRI-Mech  3.0  and  mixture-averaged  multi- 
component  transport  (details  of  computations  can  be  found  in 
[25]).  AT  represents  the  temporal  evolution  of  temperature  in 
mixture  fraction  space: 


Fig.  2  -  Influence  of  the  Soret  effects  on  the  evolution  of 
temperature  during  autoigni-tion  of  Case  H  obtained  from 
computations  of  igniting  counterflow  diffusion  flames. 
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AT(t)  =  max(T(Z,  t)  -  T(Z,  0))  (1) 

A  comparison  of  AT  between  the  computations  indicates 
that  the  contribution  of  the  Soret  effect  to  the  predicted 
chemical  time  scales  is  less  than  5%.  Such  a  small  contribu¬ 
tion  of  the  Soret  effect  in  the  ID  computations  can  be 
explained  since  this  effect  is  only  important  for  light  species 
(H,  H2)  He)  at  low  temperatures  [26].  Since  the  mass  fraction  of 
H2  in  the  fuel  stream  is  very  small  (0.0073)  and  the  tempera¬ 
ture  levels  are  high  (above  1000  K)  eventually  the  contribution 
of  the  Soret  effect  is  small  in  the  studied  conditions.  This 
observation  can  be  translated  into  the  corresponding 
computed  laminar  2D  flames  since  the  turbulence  transport  is 
absent  and  molecular  diffusion  plays  an  important  role  in  the 
flame  structure.  In  our  previous  study  [15],  it  has  been 
demonstrated  that  the  structure  of  the  2D  laminar  Mild  flames 
can  be  described  very  well  with  the  computed  ID  igniting 
flamelets.  This  is  due  to  the  autoignition  structure  of  the  Mild 
flame.  A  small  change  in  the  chemical  time  scales  of  the 
computed  ID  flames  can  hardly  influence  the  location  of  the 
maximum  heat  release  and  the  flame  structure  of  the 


and  below)  iso-contours.  First  of  all,  we  note  that  the  hydrogen 
addition  does  not  affect  the  general  flame  appearance  signif¬ 
icantly.  Indeed  the  extent  of  high  temperature  zone  of  two 
flames  is  comparable  and  the  maximum  temperatures  are 
rather  similar,  being  ~  50  K  higher  in  Case  M.  However,  the 
addition  does  appreciably  influence  the  position  of  the  tem¬ 
perature  rise  in  respect  to  the  location  of  the  burner  outlet. 
The  significant  temperature  rise  in  Case  M  is  evident  at  a 
height  of  ~6  mm,  while  that  in  Case  H  is  visible  at  a  height  of 
~  3  mm.  Such  a  behavior  indicates  that  the  hydrogen  addition 
promotes  ignition.  The  location  of  the  temperature  rise  in 
Case  H  is  also  somewhat  shifted  to  the  fuel  side  in  comparison 
with  that  in  Case  M.  We  also  note  that  the  hydrogen  addition 
decreases  significantly  the  axial  location  of  a  stoichiometric 
contour  from  ~  18  mm  to  12  mm.  The  above  comments  are 
further  illustrated  by  the  simulated  heat  release  rate  distri¬ 
butions  shown  in  Fig.  4  together  with  streamlines  and  the  iso¬ 
contours.  As  can  be  seen  from  Figs.  3  and  4,  the  location  of 
maximal  temperature  at  every  height  corresponds  to  the  po¬ 
sition  of  the  maximal  heat  release.  Fig.  4  also  shows  that  in 
both  Cases  the  combustion  starts  at  a  lean  mixture  fraction, 


computed  2D  flames.  Consequently,  the  Soret  effect  can  be 
safely  neglected  in  the  2D  computations. 


4.  Results  and  discussion 

4.1.  General  structure  of  flames 


To  illustrate  the  effects  of  hydrogen  addition  on  the  overall 
flame  structure,  we  present  two-dimensional  false  color  plots 
of  the  computed  distribution  of  temperature  for  flames  with 
and  without  hydrogen  addition  in  Fig.  3.  The  figure  includes 
the  stoichiometric  and  most  reactive  mixture  fraction  (see  [15] 


CASEM 


CASE  H 


R  (mm) 


50  - 


T  (K) 

1600 
1550 
1500 
1450 

1400  30 

1350  p" 
1300  £ 

1250  -20 
1200  ^ 
1150 
1100 


0 

R  (mm) 


T  (K) 

1550 

1500 

1450 

1400 

1350 

1300 

1250 

1200 

1150 

1100 


Fig.  3  -  Two-dimensional  false  color  plots  of  the  line 
computed  temperature  in  Cases  M  and  H.  The  black  lines 
show  the  stoichiometric,  ?st,  and  most  reactive  mixture 


proceeds  further  downstream  towards  and  beyond  the  stoi¬ 
chiometric  mixture  fraction.  We  further  note  that  the  both 
heat  release  profiles  are  notably  different  from  those  of  coflow 
diffusion  edge  flames,  see  for  examples  [15].  This  is  due  to  the 
fact  that  the  flame  temperature  in  our  Cases  is  not  large 
enough  to  support  an  edge  flame  that  can  propagate  up¬ 
stream.  Therefore,  it  can  be  concluded  that  stabilization 
mechanism  of  these  cases  is  governed  by  autoignition. 

The  autoignition  structure  of  Cases  M  and  H  is  further 
analyzed  by  plotting  the  computed  distributions  of  the  pro¬ 
duction  rate  of  formaldehyde  in  Fig.  5.  It  is  clear  that  for  both 
flames  CH20  is  formed  upstream  of  the  maximum  heat 
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Fig.  4  -  Two-dimensional  false  color  plots  of  the  computed 
heat  release  (HR)  for  Cases  M  and  H.  The  lines  with  arrows 
represent  streamlines  of  the  velocity  field.  The  black  lines 
show  the  stoichiometric  'Cst  and  most  reactive  mixture 


fraction,  i;MR,  iso-contours.  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 


fraction  Cmr  iso-contours.  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 
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Fig.  5  -  Two-dimensional  false  color  plots  of  the  computed 
production  rate  of  formaldehyde  («j(CH20))  for  Cases  M  and 
H.  The  lines  with  arrows  represent  streamlines  of  the 
velocity  field.  The  black  lines  show  the  stoichiometric  £st 
and  most  reactive  mixture  fraction  i(MR  iso-contours.  (For 
interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this 
article.) 


release  location,  immediately  above  the  burner  tip.  The  pro¬ 
duced  CH20  is  then  consumed  downstream  at  the  maximum 
heat  release  location.  We  note  that  such  a  behavior  is  in  a  clear 
contrast  with  the  trends  that  have  been  observed  for  the  edge 
flames,  see  [15],  where  CH20  is  formed  downstream  of  the 
maximum  heat  release  and  it  is  diffused  upwards  towards  the 
maximum  heat  release  location  where  it  is  consumed.  This 
again  demonstrates  the  stabilization  of  these  flames  by  the 
autoignition  phenomena.  The  location  of  the  heat  release  and 
CH20  production  rate  can  be  discussed  in  terms  of  most 
reactive  mixture  fraction  CMr  (the  mixture  fraction  with  the 
shortest  autoignition  time)  [27],  This  value  has  been  obtained 
from  the  computations  of  a  homogeneous  reactor  for  various 
frozen  mixtures  of  fuel  and  oxidizer  in  a  ID  configuration  for 
Cases  M  and  H  (more  details  can  be  found  in  [25]  and  [15]).  It  is 
observed  that  ignition  starts  from  t(MR  similar  to  what  has  been 
reported  for  the  autoignition  structures  [25,27,28].  The  t(MR  is 
located  at  the  lean  side  due  to  the  higher  temperature  and 
availability  of  a  small  amount  of  02  at  the  coflow  side. 

Figs.  3-5  also  display  the  increased  diffusivity  of  the 
mixture  due  to  hydrogen  addition.  Indeed,  since  the  methane 
and  hydrogen  were  mixed  while  keeping  the  total  flow  rate  of 
the  fuel  constant,  the  global  equivalence  ratio  has  decreased 
by  the  hydrogen  addition.  However,  the  position  of  the  stoi¬ 
chiometric  mixture  fraction  contour  near  the  burner  outlet  is 
practically  not  affected  by  the  hydrogen  addition. 

An  important  aspect  of  Case  H  compared  to  Case  M  is  the 
flame  stability  which  is  significantly  improved  by  H2  enrich¬ 
ment.  In  fact,  auto-ignition  of  these  flames  is  initiated  at  low 
mixture  fractions  (<0.12)  close  to  coflow  stream  in  localized 


regions  of  low  scalar  dissipation  rate.  Preferential  diffusion 
effects  of  H2  leads  to  an  increased  presence  of  fuel  at  this  most 
reactive  mixture  fraction  [29]  and  eventually  an  increased 
chemical  reactivity  and  stability  of  the  flame.  This  observation 
corresponds  very  well  with  recent  measurements  of  the  tur¬ 
bulent  Jet-in-Hot  Coflow  burner  of  TU  Delft  [30],  In  these 
measurements,  they  have  added  a  range  5%-25%  hydrogen  to 
the  lifted  Mild  flames  and  they  have  observed  a  significant 
improvement  in  stability  together  with  shorter  lift-off  height. 
Also  the  works  of  the  Adelaide  group  [6]  confirm  that  the  flame 
stabilization  in  MILD  combustion  is  very  similar  for  different 
fuels  (NG,  C2H4  and  LPG)  with  H2  enrichment.  Such  an 
improvement  in  the  flame  stability  induced  by  the  presence  of 
hydrogen  suggests  a  more  flexible  range  of  possible  fuels  for 
application  in  MILD  combustors. 

To  supplement  the  computational  findings  as  to  the  effect 
of  the  H2  addition  on  the  flame  structure,  we  compare  the 
radial  profiles  of  temperature  measured  at  3  and  12  mm  above 
the  burner  outlet  for  two  cases  studied  in  Fig.  6.  To  illustrate 
the  development  of  the  flames,  we  plot  the  difference  between 
the  radial  profiles  measured  in  flames  with  those  measured 
when  only  N2  was  flowing  through  the  fuel  tube.  The  situation 
when  only  N2  is  present  in  the  fuel  tube  corresponds  to  frozen 
(no  chemical  activity)  mixing  between  coflow  and  fuel.  The 
development  in  the  profiles  of  the  temperature  difference  for 
both  flames  reflect  the  “wishbone-shaped”  radial  tempera¬ 
ture  distributions  in  an  “ordinary”  diffusion  flame  [15],  with 
the  exception  that  the  temperature  rise  in  non-diluted  coflow 
is  approximately  a  factor  of  9  higher  than  that  shown  in  Fig.  6. 
The  figure  shows  that  at  3  mm  height  the  peak  in  the  tem¬ 
perature  difference  in  Case  M  is  ~  100  K,  while  that  in  Case  H  it 
is  ~  200  K.  The  experimental  data  at  3  mm  height  show  the 
presence  of  a  noteworthy  measured  amount  of  02  near  the 
centerline  in  Case  M  [15],  while  in  Case  H  the  02  mole  fraction 


Fig.  6  -  The  differences  between  the  radial  temperature 
profiles  at  3  and  12  mm  with  those  measured  when  only  N2 
was  flowing  through  the  fuel  tube.  High-order  polynomial 
trendlines  are  added  to  accentuate  the  trends  in  the 
experimental  data. 
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Fig.  7  -  Measured  (pointes)  and  computed  (lines)  radial  profiles  of  temperature  and  major  species  in  Case  H.  Red  dash  lines 
denote  the  calculations  with  a  considerably  reduced  thickness  of  the  thermal  boundary  layer.  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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is  close  to  zero  (not  shown),  similarly  to  the  data  at  2  mm  (see 
Fig.  7).  The  analysis  of  the  data  measured  in  the  vicinity  of  the 
burner  outlet  indicated  that  the  hydrogen  addition  promotes 
the  ignition.  We  would  like  also  to  mention  that  the  radial 
profiles  of  the  temperature  difference  in  Case  M  measured  at 
several  downstream  distances,  see  Fig.  6  and  [13]  show  the 
peak  in  the  temperature  difference  on  the  level  of  200  K.  This 
value  is  similar  to  those  for  Case  H.  The  figure  also  shows  that 
the  position  of  the  maximal  heat  release  shifts  notably  to¬ 
wards  the  burner  outlet  with  the  addition.  We  would  like  to 
note  that  some  of  the  described  above  effects  of  the  hydrogen 
addition  on  general  flame  structure  such  as  reduction  of  the 
flame  height  and  the  shift  of  the  reaction  zone  towards  the 
burner  exit  were  also  reported  in  [8]. 

4.2.  Comparison  of  the  experimental  results  measured  in 

Case  H  with  the  numerical  simulations 

As  an  example  of  the  measurements  performed,  Fig.  7  shows 
the  experimental  radial  profiles  of  temperature  and  major 
species  mole  fractions  at  the  axial  distances  2,  6  and  12  mm  in 
Case  H.  To  test  the  predictive  power  of  the  GRI  3.0  mechanism 
at  the  specific  experimental  conditions,  we  contrast  the 
measurements  with  the  results  of  the  computation.  The 
progress  of  combustion  in  Case  H  is  reflected  in  the  develop¬ 
ment  of  temperature.  However  the  details  of  the  evolution  is 
considerably  masked  by  the  large  temperature  difference  be¬ 
tween  the  coflow  and  fuel  similarly  as  in  the  flame  without 
hydrogen  addition  (Case  M),  see  [13,15],  The  development  of 
the  major  species  profiles  for  Case  H  is  consistent  with  tem¬ 
perature  changes.  The  experimental  temperature  agrees 
quantitatively  with  the  computations  at  the  axial  distance  of 
2  mm.  We  note  that  the  little  bump  in  the  experimental 
temperature  profile  (marking  apparently  the  location  of 
combustion  zone)  at  the  radial  distance  of  ~5  mm  is  well 
captured.  At  the  axial  distances  of  6  and  12  mm,  the  agree¬ 
ment  is  also  good,  the  calculations  predict  the  measured  data 
within  120  K  and  generally  better.  The  calculations  predict 
rather  well  the  H2  profiles  and  somewhat  overpredict  the  CH4 
consumption.  The  02  profile  is  quantitatively  predicted  at  the 
height  of  2  mm.  We  noted  already  the  presence  of  significant 
amount  of  02  near  the  centerline.  Further  downstream, 
however,  the  decrease  in  the  experimental  02  fraction  at  the 
edges  of  the  measured  domain  occurs  notably  closer 
(~  1.5  mm)  to  the  coflow  side  than  that  in  calculations.  Such  a 
behavior  implies  that  the  computational  combustion  zone  is 
shifted  in  comparison  with  the  experimental  one  closer  to  the 
fuel  side.  The  overprediction  of  C02  and  H20  profiles  inside  the 
domain  of  the  fuel  flow  for  some  axial  distances  is  consistent 
with  this  shift. 

Fig.  8  shows  the  experimental  and  calculated  vertical  pro¬ 
files  of  T  and  CH4  and  H2  mole  fractions.  The  experimental 
temperature  rapidly  increases  to  the  maximum  temperature 
of  ~1580  K,  at  the  distance  of  ~15  mm  above  the  burner 
outlet.  The  computations  support  the  maximum  measured 
temperature  quantitatively,  however,  approach  it  faster  and 
reach  the  maximum  already  at  a  distance  of  ~12  mm.  In 
accordance  with  the  observed  difference  in  the  behavior  of  the 
temperature  profiles,  the  calculations  and  measurements 
show  similar  discrepancy  in  the  fuel  profiles.  While  the 
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Fig.  8  -  The  experimental  and  calculated  vertical  profiles  of 
temperature  and  CH4  and  H2  mole  fractions  in  Case  H. 


experimental  data  demonstrate  the  non-negligible  CH4  and  H2 
mole  fractions  at  distances  above  12  mm,  the  computations 
show  that  all  fuel  is  consumed  by  a  height  of  12  mm.  The  rapid 
initial  decrease  in  the  vertical  CH4  and  H2  profiles  is  a  result  of 
the  diffusion  of  fuel  in  radial  directions  towards  the  flame 
front.  We  would  like  to  point  out  to  the  interesting  feature, 
qualitatively  well  captured  also  by  the  calculations,  in 
behavior  of  the  vertical  CH4  and  H2  profiles  in  the  vicinity  of 
the  maximal  heat  release.  The  CH4  concentration  is  consumed 
rapidly  in  the  fuel-rich  part  of  the  flame  front  (see  also  Fig.  9), 
while  most  of  the  H2  concentration  survives  this  region;  the 
hydrogen  burns  out  completely  further  downstream  at  the 
stoichiometric  and  fuel  lean  conditions.  The  behavior  is  not 
perhaps  unexpected,  the  one-dimensional  flame  calculations 
of  the  CH4/H2  blends  demonstrate  that  a  significant  part  of  the 


Fig.  9  -  The  vertical  profiles  of  mixture  fraction  and  heat 
release  (HR)  in  Case  H  (symbols:  measurements,  lines: 
computations).  The  horizontal  line  marks  the  position  of 
stoichiometric  mixture  fraction. 
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hydrogen  is  not  consumed  at  fuel-rich  conditions.  Fig.  9  shows 
the  vertical  profiles  of  heat  release  and  the  mixture  fractions. 
The  mixture  fraction  was  calculated  using  Bilger’s  formula 
[31]  modified  to  account  for  the  presence  of  C  and  H  in  the 
oxidizer.  The  horizontal  line  in  Fig.  9  marks  the  position  of 
stoichiometric  mixture  fraction.  Here  we  point  out  that  the 
mixture  fraction  is  predicted  well,  but  the  CH4  and  H2  profiles 
are  less  accurately  reproduced  (Fig.  8). 

Although  the  computations  in  Case  H  reproduce  quite  well 
the  general  trend  of  the  experimental  data,  the  results  of  this 
section  indicated  some  quantitative  differences.  The  exami¬ 
nation  of  the  behavior  of  the  experimental  and  calculated 
profiles  suggests  that  the  computational  flame  burns  faster 
than  the  experimental  one. 

A  possible  reason  for  the  differences  might  be  the  choice  of 
the  boundary  conditions.  Namely,  in  Cases  M  and  H,  the 
thermal  boundary  layer  between  the  fuel  tube  and  the  coflow 
was  chosen  in  a  way  to  reproduce  the  measured  values  of 
temperature  for  the  non-reacting  mixture  [15].  To  do  this,  CH4 
has  been  completely  removed  and  replaced  by  N2,  keeping  all 
other  parameters  of  the  Case  M  the  same,  to  avoid  any 
possible  heat  release.  The  calculated  temperature  corre¬ 
sponded  very  well  with  the  measured  values  at  a  height  of 
3  mm  and  12  mm  above  the  fuel  tube  exit  and  further 
modeling  of  Case  M  provided  a  very  good  match  with  the 
experimental  data.  In  [15],  however,  we  did  not  check  the 
sensitivity  of  the  numerical  results  to  the  chosen  boundary 
conditions.  In  an  attempt  to  analyze  the  sensitivity,  we  also 
performed  the  numerical  calculations  with  a  considerably 
reduced  thickness  of  the  thermal  boundary  layer  (see  [15]). 
Fig.  7  includes  the  radial  profiles  of  temperature  and  species 
concentrations  calculated  using  these  new  boundary  condi¬ 
tions  at  a  height  of  6  mm.  The  figure  shows  that  decreasing  the 
boundary  layer  thickness  does  not  affect  significantly  both  the 
temperature  and  species  distributions. 

While  we  can  not  exclude  other  possible  reasons  for  the 
differences  observed  (for  example,  uncertainties  in  kinetics 
and  transport,  or  uncertainty  in  vertical  position),  we  believe 
that  a  further  analysis  of  the  physical/chemical  reasons  would 
hardly  be  justified  and  is  beyond  the  scope  of  the  manuscript. 
We  are  also  of  the  opinion  that  such  an  analysis  would  require 
an  addition  experimental  efforts,  perhaps,  the  measurements 
of  the  OH  radical  distribution  might  be  beneficial  for  this 
purpose  (we  plan  to  conduct  such  measurements  in  the 


Fig.  11  -  Measured  and  calculated  NO  mole  fractions 
versus  mixture  fraction  in  Case  H.  The  vertical  line  marks 
the  position  of  stoichiometric  mixture  fraction. 


future).  Furthermore,  there  may  not  exist  a  single  reason 
which  explains  these  differences. 

4.3.  The  NO  formation  in  CH4/H2  flame 

Fig.  10  displays  the  NO  radial  profiles  measured  and  calculated 
at  axial  positions  of  3,  6  and  15  mm  in  Case  H.  At  the  axial  dis¬ 
tances  of  3  and  6  mm  the  NO  profiles  show  a  minimum  at  the 
flame  centre,  at  15  mm  the  profile  becomes  flat.  Comparison 
with  the  calculated  profiles  shows  excellent  agreement.  We 
remind  here  that  a  similar  trend  in  the  behavior  of  the  NO  radial 
profile  with  downstream  distance  was  noticed  in  Case  M  [13,14], 
In  these  papers  it  was  also  found  (by  analyzing  the  experimental 
and  computational  data  in  mixture  fraction  space)  that  the  NO 
field  is  dominated  by  mixing  of  the  NO  formed  in  the  coflow  with 
the  reaction  products  of  the  diluted  fuel,  with  negligible  NO 
formation  from  the  fuel  via  the  Fenimore  mechanism.  To  check 
the  effects  of  H2  addition  on  the  NO  formation,  we  plot  the 
experimental  and  computational  results  as  a  function  of 


Z  =  3  mm  Z  =  6  mm  Z  =  1 5  mm 


Fig.  10  —  The  measured  and  calculated  radial  profiles  of  NO  in  Case  H  at  axial  positions  of  3,  6  and  15  mm.  Pointes  denote 
measurements;  lines  denote  calculations. 
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mixture  fraction  in  Case  H,  see  Fig.  11.  The  figure  also  includes 
the  mixing  curve  indicating  simple  mixing  of  the  hot  coflow 
with  the  combustion  products  of  the  fuel  jet,  (NO)  =  (1  -  £) 
(NO)oxWm/Wox,  where  £  is  the  mixture  fraction,  (NO)ox  is  NO 
fraction  in  the  oxidizer  (10  ppm)  and  WM  and  Wox  are  molecular 
weights  of  the  combustion  products  at  the  measuring  position 
and  oxidizer  stream,  respectively.  The  measurements  and  cal¬ 
culations  at  3  and  15  mm  cluster  around  the  mixing  curve 
demonstrating  almost  no  NO  formation  or  consumption.  The 
comparison  of  the  measured  NO  fraction  at  6  mm  with  the 
mixing  line  demonstrates  a  noticeable  NO  consumption  at  the 
fuel  rich  conditions,  the  difference  is  on  the  level  ~  30%,  clearly 
outside  the  measurement  uncertainty  of  20%.  We  note  that  in 
Case  M  at  the  given  level  of  the  NO  seeding  in  the  coflow,  the 
difference  was  on  the  level  of  20%.  At  the  fuel  lean  side  the  date 
concentrate  around  the  mixing  line.  The  calculations  support 
the  measurements  quantitatively,  although  they  predict  a 
somewhat  smaller  difference  with  the  mixing  line  ( —  20%)  at 
fuel  rich  conditions.  The  rate-of-production  analysis  reveals 
that  the  consumption  is  mostly  due  to  the  conversion  of  the  NO 
molecules  to  the  N-containing  intermediates  (HCN,  NH3  and 
HNCO).  These  intermediates  (similarly  to  Case  M)  are  converted 
back  to  NO  at  downstream  distances. 


5.  Conclusions 

An  experimental  and  computational  study  of  the  effects  of  the 
H2  addition  on  both  the  spatial  structure  and  the  NO  formation 
of  a  MILD  flame  burning  in  an  axisymmetric  laminar-jet-in- 
hot-coflow  geometry  has  been  performed.  Using  methane  as 
base  fuel,  the  effects  of  adding  H2  were  studied.  The  addition  of 
hydrogen  appreciably  decreases  the  flame  height  ( — 25%)  and 
promotes  the  ignition,  however  only  modestly  affects  the 
maximal  flame  temperature  and  the  thickness  of  combustion 
zone.  The  numerical  data  are  in  reasonably  good  agreement 
with  the  measurements.  Analysis  of  the  NO  fraction  as  a 
function  of  mixture  fraction  indicates  that  the  NO  distribution 
is  dominated  by  mixing  of  the  NO  formed  in  the  coflow  with  the 
reaction  products  of  the  diluted  fuel,  with  negligible  NO  for¬ 
mation  from  the  fuel  in  all  flames  studied.  Experiments  show 
the  considerable  reduction  of  the  seeded  NO  at  a  height  of 
6  mm.  The  calculations  support  the  observation  quantitatively 
and  reveal  that  the  reduction  is  due  to  the  conversion  of  the  NO 
molecules  to  the  N-containing  intermediates  which  are  con¬ 
verted  back  to  NO  at  downstream  distances. 
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